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Continental  runoff  is  a  major  source  of  freshwater,  nutrients  and  terrigenous  material  to  the  Arctic  Ocean. 
As  such,  it  influences  water  column  stratification,  light  attenuation,  surface  heating,  gas  exchange,  biological 
productivity  and  carbon  sequestration.  Increasing  river  discharge  and  thawing  permafrost  suggest  that  the 
impacts  of  continental  runoff  on  these  processes  are  changing.  Here,  a  new  optical  proxy  was  developed  and 
implemented  with  remote  sensing  to  determine  the  first  pan- Arctic  distribution  of  terrigenous  dissolved 
organic  matter  (tDOM)  and  continental  runoff  in  the  surface  Arctic  Ocean.  Retrospective  analyses  revealed 
connections  between  the  routing  of  North  American  runoff  and  the  recent  freshening  of  the  Canada  Basin, 
and  indicated  a  correspondence  between  climate-driven  changes  in  river  discharge  and  tDOM  inventories  in 
the  Kara  Sea.  By  facilitating  the  real-time,  synoptic  monitoring  of  tDOM  and  freshwater  runoff  in  surface 
polar  waters,  this  novel  approach  will  help  understand  the  manifestations  of  climate  change  in  this  remote 
region. 


Continental  runoff  influences  surface  ocean  processes  throughout  the  Arctic.  Rivers  discharge  about 
3,300  km3  yr_1  of  buoyant  freshwater  to  the  Arctic  Ocean1  where  it  sustains  a  halocline  that  isolates 
the  upper  200-m  layer  from  a  warmer,  saltier  Atlantic  layer2.  This  upper  layer  represents  only  0.1%  of  the 
global  ocean  volume  but  receives  11%  of  global  riverine  discharge  from  a  circumpolar  drainage  basin  that 
encompasses  over  10%  of  the  world’s  land  surface  area.  Riverine  input  is  evident  throughout  the  Arctic  from 
low  salinities  and  high  concentrations  of  terrigenous  dissolved  organic  matter  (tDOM)  in  polar  surface  waters3,4. 
Arctic  rivers  discharge  onto  the  world’s  largest  continental  shelves  that  surround  and  cover  more  than  half  of  the 
Arctic  Ocean  area5.  These  river-influenced  ocean  margins  are  dominant  features  of  the  Arctic  and  important 
regions  of  sea-ice  production,  biological  productivity,  carbon  sequestration  and  processes  that  drive  biogeochem¬ 
ical  cycles  on  a  basin-wide  scale. 

The  impact  of  rivers  on  the  freshwater  content  and  biogeochemistry  of  the  Arctic  Ocean  is  strongly  influenced 
by  the  residence  time  of  riverine  inputs  and  their  pathways  to  the  North  Atlantic  Ocean.  The  average  residence 
time  of  Arctic  shelf  water  ranges  from  months  to  years  depending  on  shelf  region6,7.  Variations  in  atmospheric 
forcing  can  affect  residence  time  by  dramatically  altering  continental  runoff  routing  and  distribution  in  the  Arctic 
Ocean7,8.  Climate-driven  changes  in  runoff  distribution  have  been  linked  to  large  inter-annual  and  decadal 
variations  in  the  freshwater  content  of  the  Arctic  Ocean,  and  these  fluctuations  have  important  ramifications 
for  North  Atlantic  meridional  overturning  circulation9.  Furthermore,  changes  in  runoff  routing  through  the 
Arctic  Ocean  can  regulate  the  extent  to  which  Arctic  tDOM  is  incorporated  into  North  Atlantic  Deep  Water 
(NADW)  and  distributed  in  the  global  ocean4.  Despite  their  global  significance,  continental  runoff  routing  and 
distribution  in  the  Arctic  remain  poorly  understood  and  difficult  to  predict  in  this  remote  and  under-sampled 
region  of  the  global  ocean.  Here,  a  new  optical  proxy  for  tDOM  implemented  using  remote  sensing  of  ocean  color 
facilitates  the  synoptic  and  retrospective  monitoring  of  continental  runoff  in  the  Arctic  Ocean,  and  provides  a 
novel  approach  to  assess  climate-driven  changes  in  the  Arctic. 
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Results 

Dissolved  lignin  is  well  established  as  a  biomarker  of  tDOM  in 
oceanic  waters10,11  and  has  been  successfully  applied  as  a  tracer  of 
riverine  inputs  in  the  Arctic  Ocean3,4.  However,  its  use  relies  on  the 
collection  of  seawater  samples  and  land-based  laboratory  analyses 
that  prevent  synoptic  applications.  Lignin  is  also  an  important  chro- 
mophore  in  tDOM,  a  property  that  facilitates  detection  using  optical 
properties.  Here,  we  demonstrate  that  the  spectral  slope  coefficient  of 
chromophoric  dissolved  organic  matter  (CDOM)  between  2 75  and 
295  nm  (S275-295)  is  a  practical  and  reliable  optical  proxy  for  dis¬ 
solved  lignin  and  tDOM  across  various  river-influenced  ocean  mar¬ 
gins  of  the  Arctic  Ocean12,13  (Fig.  la).  An  increase  in  S275-295  is 
indicative  of  a  decrease  in  the  tDOM  content  of  surface  waters. 
Under  the  assumption  that  S275-295  behaves  conservatively,  an 
increase  in  $275-295  corresponds  to  a  diminishing  influence  of  con¬ 
tinental  runoff.  Biological  degradation  has  little  effect  on  $275-295 13,14 
and  photodegradation  of  tDOM  in  polar  surface  waters  is  not  extens¬ 
ive4,15,  thereby  warranting  the  assumption  and  the  use  of  $275-295  as  a 
tracer  of  tDOM  and  freshwater  runoff  in  river- influenced  ocean 
margins  of  the  Arctic. 

The  practicality  of  the  $275-295  proxy  resides  in  its  amenability  to 
remote  sensing.  The  development  of  an  empirical  algorithm  reveals 
that  $275-295  can  be  estimated  from  ocean  color  with  excellent  and 
consistent  accuracy  (within  ±  4%)  over  the  broad  range  of  water 
types  and  $275-295  values  observed  in  the  Arctic  Ocean  (Fig.  lb).  The 
strong  connection  between  $275-295  and  ocean  color  in  river-influ¬ 
enced  ocean  margins  exists  because  the  optical  property  $275-295  is  an 
excellent  tracer  of  tDOM13  and  the  gradients  of  ocean  color  in  river- 
influenced  ocean  margins  are  primarily  driven  by  the  constituents  of 
continental  runoff.  The  $275-295  algorithm  thus  facilitates  the  syn¬ 
optic  and  real-time  monitoring  of  tDOM  and  freshwater  runoff  in 
the  Arctic  using  ocean-color  remote  sensing. 

Implementation  of  the  $275-295  algorithm  using  MODIS  Aqua 
ocean  color  provides  the  first  pan-Arctic  view  of  average  tDOM 


distributions  in  river- influenced  ocean  margins  during  August 
2002-2009  (Fig.  lc).  The  influence  of  the  three  largest  Arctic  rivers 
on  tDOM  distributions  is  evident  along  the  entire  Siberian  margin. 
One  of  the  most  striking  features  is  the  overwhelming  amount  of 
tDOM  in  the  Laptev  and  East  Siberian  Seas,  indicating  a  large  portion 
of  Eurasian  river  tDOM  and  runoff  is  routed  alongshore  in  an  east¬ 
ward  direction.  This  observation  is  supported  by  measurements  of 
high  in  situ  concentrations  of  dissolved  organic  carbon  of  terrigenous 
origin  in  the  East  Siberian  Sea  in  August  2008 16,  and  is  consistent 
with  shelf  water  residence  times  of  several  years  in  this  region6,7.  Low 
$275-295  values  along  the  Alaskan  coast  in  the  Chukchi  Sea  depict  the 
transport  of  Yukon  River  tDOM  to  the  Arctic  Ocean  in  the  Alaska 
Coastal  Current  and  indicate  notable  tDOM  contributions  by  smal¬ 
ler  Alaskan  rivers  (e.g.,  Noatak  and  Kobuk  Rivers).  Most  of  the  runoff 
and  tDOM  from  the  Mackenzie  River  enters  the  Beaufort  Sea  and 
appears  to  be  entrained  in  the  gyre  circulation  of  the  southern 
Canada  Basin.  A  relatively  small  fraction  of  the  discharge  from  the 
Mackenzie  River  was  routed  to  the  Labrador  Sea  via  the  Canadian 
Arctic  Archipelago  (CAA). 

Assessment  of  inter- annual  variability  in  tDOM  distribution 
revealed  a  change  in  the  routing  of  Mackenzie  River  discharge  over 
the  last  decade  (Fig.  2;  Supplementary  Figs.  S3  and  S4).  The 
Mackenzie  River  outflow  progressed  from  an  alongshore,  eastward 
path  to  the  CAA  in  2002  (Fig.  2a)  to  a  cross-shelf,  northwestward 
path  to  the  Canada  Basin  since  2006  (Fig.  2b).  This  routing  switch  has 
important  implications  for  the  fate  of  North  American  runoff  and 
tDOM.  Routing  through  the  CAA  favors  the  export  of  North 
American  tDOM  to  the  Labrador  Sea17.  Although  the  magnitude  of 
this  export  remains  unknown,  deep  ocean  convection  in  the 
Labrador  Sea  can  sequester  tDOM  in  North  Atlantic  Deep  Water 
for  centuries4.  In  contrast,  routing  to  the  Canada  Basin  contributes  to 
longer-term  storage  of  runoff  and  processing  of  tDOM  within  the 
Beaufort  Gyre  and  Arctic  Ocean4,15.  This  decadal  shift  to  greater 
export  of  Mackenzie  River  runoff  to  the  Canada  Basin  during 
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Figure  1  |  The  CDOM  spectral  slope  coefficient,  $275-295?  as  a  tracer  of  riverine  inputs  in  the  Arctic  Ocean,  (a)  The  relationship  between 
$275-295  and  dissolved  lignin  carbon  yield  (TDLP9-C)  across  various  Arctic  river-influenced  ocean  margins,  (b)  A  performance  evaluation  of  the  $275-295 
algorithm.  On  average,  $275-295  can  be  estimated  from  ocean  color  within  4%  of  $275-295  values  measured  in  situ,  (c)  Implementation  of  the 
algorithm  using  MODIS  Aqua  ocean  color  providing  a  pan-Arctic  view  of  an  August  climatology  (2002-2009)  of  $275-295-  Increasing  $275-295  values  are 
indicative  of  a  decreasing  fraction  of  tDOM.  The  four  largest  Arctic  rivers  are  labeled  and  ranked  in  order  of  decreasing  discharge:  Yenisei  (1),  Lena  (2), 
Ob  (3),  and  Mackenzie  (4).  River-influenced  margins  of  the  Arctic  are  labeled:  Gulf  of  Ob  (GO),  Kara  Sea  (KS),  Laptev  Sea  (LS),  East  Siberian  Sea 
(ESS),  Chukchi  Sea  (CS),  Beaufort  Sea  (BS)  and  Amundsen  Gulf  (AG).  The  contour  lines  represent  the  2000-m  isobath. 
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Figure  2  |  Inter-annual  changes  in  the  routing  of  Mackenzie  River  runoff,  (a)  The  Mackenzie  River  outflow  in  2002,  showing  eastward  transport  of 
North  American  runoff  to  the  Canadian  Arctic  Archipelago,  (b)  The  Mackenzie  River  outflow  in  2008,  showing  northwestward  transport  of  North 
American  runoff  to  the  Canada  Basin.  A  progressive  switch  from  eastward  to  northwestward  routing  occurred  between  2002  and  2011  (Supplementary 
Figs.  S3  and  S4)  and  coincides  with  the  rapid  freshening  of  the  Canada  Basin.  The  S2 75-295  algorithm  was  implemented  using  4-km  resolution,  yearly- 
binned,  MODIS  Aqua  ocean  color.  The  Mackenzie  River  is  labeled  (4).  The  contour  line  represents  the  2000-m  isobath  and  outlines  the  Canada  Basin. 


2006-2011  could  contribute  to  less  North  American  tDOM  being 
sequestered  in  the  deep  ocean. 

The  routing  of  Mackenzie  River  outflow  since  2006  coincides  with 
observations  of  a  rapid  accumulation  of  freshwater  in  the  Canada 
Basin.  Freshening  began  in  the  1990s  and  accelerated  in  the  late 
2000s18-21.  The  realization  that  a  large  accumulation  of  freshwater 
in  the  Canada  Basin  could  impact  global  ocean  circulation  stimulated 
research  to  identify  the  freshwater  sources  and  climatic  forcing 
responsible  for  observed  changes  in  salinity.  Pacific  water,  ice-melt, 
precipitation  and  river  runoff  are  distinct  sources  of  freshwater  to  the 
Canada  Basin.  The  relative  contributions  of  these  freshwater  sources 
to  the  freshening  are  typically  estimated  using  tracers  (e.g.  dlsO, 
alkalinity,  salinity)  or  are  indirectly  calculated  using  complex 


models9,20,22-25.  The  lignin-based  $275-295  Proxy  presented  here  sub¬ 
stantiates  the  riverine  origin  of  freshwater  and  provides  a  direct 
means  to  distinguish  continental  runoff  from  other  freshwater 
sources  using  remote  sensing. 

The  acceleration  of  the  Beaufort  Gyre  due  to  the  strengthening  of 
the  Beaufort  High26  and  routing  of  Eurasian  river  runoff  to  the 
Canada  Basin  driven  by  the  Arctic  Oscillation9  are  complementary 
and  currently  leading  explanations  for  the  freshening  of  the  Basin27. 
The  remotely  sensed  runoff  distributions  presented  in  this  study 
(Fig.  2  and  Supplementary  Figs.  S3  and  S4)  further  suggest  the  recur¬ 
rent  influx  of  Mackenzie  River  runoff  to  the  Canada  Basin  in  the  late 
2000s  was  also  a  significant  source  of  freshwater  to  the  Basin  and 
contributed  to  the  freshening.  The  temporal  sequence  of  the 
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Figure  3  |  Temporal  sequence  of  the  Mackenzie  River  freshet  release  into  the  Canada  Basin  during  June  and  July  2008.  The  $275-295  algorithm  was 
implemented  using  4-km  resolution,  8-day-binned,  MODIS  Aqua  ocean  color  acquired  between  June  1  and  August  3.  The  discharge  value  (km3) 
displayed  on  each  panel  indicates  the  cumulative  Mackenzie  River  discharge  released  between  May  1  (start  of  the  freshet)  and  the  end  of  each  8-day  period. 
The  contour  line  represents  the  2000-m  isobath  and  delineates  the  Canada  Basin.  The  Mackenzie  River  is  labeled  (4). 
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Figure  4  |  Climatic  forcing  of  major  Eurasian  riverine  inputs,  (a)  The  North  Atlantic  Oscillation  (NAO)  and  Arctic  Oscillation  (AO)  wintertime  indices 
(DJFM:  Dec-Jan-Feb-Mar)  in  relation  to  the  combined  discharge  anomaly  of  the  Yenisei,  Lena  and  Ob  rivers  between  2002  and  2009.  High  NAO  and  AO 
indices  are  typically  associated  with  higher  precipitation  over  Northern  Eurasia.  River  discharge  is  calculated  as  the  cumulative  discharge  from  May  1  to 
August  30,  and  the  anomaly  is  relative  to  the  2002-2009  average,  (b)  Average  S27 5-295  values  for  surface  waters  of  the  Kara  Sea  in  relation  to  the  combined 
discharge  of  the  Yenisei  and  Ob  rivers  between  2002  and  2009.  The  $275-295  values  for  the  Kara  Sea  were  calculated  as  the  average  value  in  the  region 
displayed  in  panel  (c).  River  discharge  is  calculated  as  the  cumulative  discharge  from  May  1  to  August  30.  (c)  Distributions  of  riverine  inputs  during  high 
(2002)  and  low  (2005)  discharge  years.  The  years  correspond  to  low  and  high  average  $27 5_295  values  observed  in  the  Kara  Sea  between  2002  and  2009.  The 
S2 75_295  algorithm  was  implemented  using  MODIS  Aqua  ocean  color  data  averaged  from  May  to  October.  The  Yenisei  and  Ob  rivers  are  labeled  (1)  and 
(3),  respectively. 


Mackenzie  River  discharge  in  2008  provides  unambiguous  evidence 
that  the  spring  freshet  was  directly  and  rapidly  injected  into  the 
Canada  Basin  during  June  and  July  (Fig.  3).  In  situ  measurements 
of  salinity,  dlsO,  and  alkalinity  from  an  independent  study20  corrob¬ 
orate  that  Mackenzie  River  runoff  was  a  minor  source  of  freshwater 
to  the  Canada  Basin  in  the  early  2000s  but  increased  in  the  southern 
part  of  the  Basin  in  2007. 

The  contribution  of  North  American  runoff  to  the  Canada  Basin 
freshwater  budget  appears  to  be  shifting  in  response  to  climatic  for¬ 
cing.  A  freshwater  budget  based  on  data  from  2003-2004  estimated 
about  40%  of  North  American  runoff  enters  the  Canada  Basin22.  The 
mean  freshwater  residence  time  in  the  Basin  is  about  10  years22,  so 
this  estimate  is  representative  of  the  decade  leading  up  to  2004. 
Although  significant  freshwater  contributions  from  the  Mackenzie 
River  to  the  Basin  were  observed  during  that  decade28,  this  estimate 
does  not  integrate  the  increasing  contributions  from  North  Ameri¬ 
can  runoff  since  2006.  The  detailed  progression  of  the  2008 
Mackenzie  River  outflow  indicated  up  to  155  km3  of  Mackenzie 
freshwater  was  released  into  the  Canada  Basin  in  June  and  July  alone 
(Fig.  3),  indicating  more  than  half  of  the  discharge  from  the  Mac¬ 
kenzie  has  been  entering  the  Canada  Basin  since  2006.  Increasing 
contributions  of  North  American  runoff  to  the  Canada  Basin  are 


consistent  with  enhanced  Ekman  pumping  resulting  from  the 
wind-driven  spin-up  of  the  Beaufort  Gyre26  and  therefore  appear 
to  result  from  recent  trends  in  climate  variability. 

Discharge  from  major  Arctic  rivers  has  increased  over  the  past 
decades  in  response  to  climate  variability29-31.  Climatic  forcing  of 
river  discharge  is  evident  between  2002  and  2009  from  the  strong 
correspondence  between  the  North  Atlantic  Oscillation  (NAO)  or 
Arctic  Oscillation  (AO)  wintertime  indices  and  the  river  discharge 
anomaly  of  Eurasian  rivers  (Fig.  4a).  Manifestations  of  inter¬ 
annual  discharge  fluctuations  in  the  Arctic  Ocean  are  best  illu¬ 
strated  in  the  Kara  Sea,  which  receives  inputs  from  the  Ob  and 
Yenisei  rivers.  The  Kara  Sea  is  a  semi-enclosed  margin  that  is 
largely  ice-free  during  summer,  which  facilitates  inter-annual 
comparisons  of  remotely  sensed  variables.  Inter-annual  compar¬ 
isons  revealed  consistency  between  discharge  from  the  Ob  and 
Yenisei  rivers  and  the  average  $275-295  value  f°r  the  Kara-Sea 
sector  (Fig.  4b).  The  correspondence  between  discharge  and 
tDOM  in  the  Kara  Sea  is  particularly  striking  when  comparing 
2002  and  2005,  which  are  representative  of  high  and  low  discharge 
years  between  2002  and  2009  (Fig.  4c).  This  demonstrates  climate- 
driven  changes  in  the  tDOM  inventory  of  the  surface  Arctic 
Ocean  are  traceable  using  this  approach. 


SCIENTIFIC  REPORTS  |  3  :  1053  |  DOI:  1 0.1 038/srep01 053 


4 


Discussion 

The  S275-295  proxy  and  algorithm  presented  in  this  study  facilitate  the 
real-time  and  synoptic  monitoring  of  tDOM  and  freshwater  runoff 
in  surface  polar  waters.  The  foundation  of  the  $275-295  proxy  on  the 
lignin  biomarker  is  fundamental  to  the  approach.  The  “terrigenous” 
attribute  of  lignin  substantiates  the  riverine  origin  of  dissolved 
organic  matter  and  freshwater.  Its  “chromophoric”  property,  on 
the  other  hand,  facilitates  implementation  with  remote  sensing. 
Although  the  Arctic  remains  one  of  the  most  challenging  environ¬ 
ments  for  remote  sensing,  it  is  also  a  region  where  remote  sensing  has 
much  to  offer  given  the  logistical  restrictions  and  high  costs  of  field 
surveys,  particularly  on  the  Siberian  shelves.  The  applicability  of  the 
approach  will  improve  as  the  extent  and  duration  of  sea-ice  cover 
decrease,  allowing  for  greater  spatial  and  temporal  coverage. 

New  capabilities  to  monitor  surface  tDOM  distributions  will  prove 
valuable  for  understanding  future  climate-driven  changes  in  the  bio¬ 
geochemistry  of  the  Arctic.  Permafrost  thawing  and  precipitation  in 
drainage  basins  of  Arctic  rivers  are  projected  to  intensify  with  escal¬ 
ating  atmospheric  temperatures  and  thereby  enhance  the  mobiliza¬ 
tion  of  soil  organic  matter  to  the  Arctic  Ocean32-36.  The  fate  of  this 
material  and  its  effects  on  biogeochemical  cycles  will  depend  on  its 
routing  and  residence  time  in  polar  surface  waters4.  Furthermore, 
processing  of  organic  matter  in  surface  waters  is  influenced  by  factors 
that  are  currently  being  altered  by  climate  change  (e.g.,  ice  cover, 
water  temperature).  Multiyear  records  of  remotely  sensed  tDOM 
distributions  provide  direct  evidence  of  how  the  routing,  inventory, 
storage  and  residence  time  of  tDOM  in  surface  polar  waters  change 
in  response  to  climatic  forcing.  Such  applications  will  thus  provide 
crucial  information  for  understanding  the  fate  of  immobilized  soil 
organic  matter  and  its  impacts  on  biogeochemical  cycles  in  this 
rapidly  changing  Arctic  environment. 

Finally,  remote  sensing  of  continental  runoff  provides  useful 
insights  about  the  sources  of  freshwater  in  the  Arctic  Ocean.  Ice  melt, 
precipitation  and  runoff  are  all  increasing  under  the  current  climatic 
trends31,37,  and  are  altering  the  freshwater  budget  of  the  Arctic  Ocean. 
Remotely  sensed  runoff  distributions  provide  direct  evidence  of  con¬ 
tinental  runoff  contributions  to  specific  regions  of  the  surface  Arctic 
Ocean,  and  will  thereby  help  understand  the  mechanisms  respons¬ 
ible  for  future  changes  in  the  Arctic  freshwater  budget.  It  is  important 
to  remember,  however,  that  the  uncoupling  of  river  water  from  its 
tDOM  component  during  the  winter  freeze-thaw  cycle38  hampers  the 
use  of  this  approach  to  assess  year-to-year  changes  in  freshwater 
runoff  storage.  Future  applications  will  also  likely  include  studying 
the  influence  of  continental  runoff  on  biological  processes,  such  as 
primary  production,  in  polar  surface  waters. 

Methods 

Surface  water  samples  and  data  ( n  =  106)  used  to  establish  the  relationship 
between  S275-295  and  TDLP9-C  (Fig.  la)  were  acquired  during  several  field 
expeditions  to  the  Arctic  Ocean  between  2005  and  2010,  including  the  Amundsen 
Gulf,  Gulf  of  Ob,  and  the  Beaufort,  Chukchi,  Laptev  and  East  Siberian  Seas 
(Supplementary  Fig.  SI).  Samples  for  DOC  and  CDOM  were  filtered  onboard  and 
analyzed  in  the  home  laboratory  using  high  temperature  combustion39  and  a  dual¬ 
beam  UV-visible  spectrophotometer12,13,  respectively.  Samples  for  dissolved  lignin 
(1-10  L)  were  extracted  onboard40  and  analyzed  using  the  cupric  oxide  oxidation 
method  and  gas  chromatography  with  mass  spectrometry41.  DOC-normalized 
yields  of  nine  lignin  phenols  (TDLP9-C)  were  calculated13,  and  the  non-linear 
regression  of  TDLP9-C  on  S27 5-295  utilized  the  exponential  function  shown  in 
equation  (1), 

TDLP9— C=  exp(fl  + ^275—295)  +  exp(yS275— 295)  +  ^exp(S275— 295)  (1) 

where  a  =  2.7428,  /3  =  —245.5899,  y  =  1.6669,  and  <5  =  —0.9932  are  the  derived 
regression  coefficients. 

The  in  situ  measurements  of  ocean  color  ( n  =  236)  and  concurrent  CDOM 
samples  used  in  the  S275-295  algorithm  development  (Fig.  lb)  were  collected  during 
the  Malina  and  ICESCAPE  field  expeditions  (Supplementary  Fig.  S2).  Multispectral 
remote-sensing  reflectances,  Rrs(/ 1),  were  derived  for  5  wavelengths  (2  =  443,  488, 
531,  555  and  667  nm)  from  profiles  made  with  a  Biospherical  Compact- Optical 
Profiling  System  (C-OPS)42,  and  the  corresponding  CDOM  samples  were  collected 
and  analyzed  according  to  NASA  standard  protocols43.  “Measured”  S275-295  were 


calculated  from  CDOM  absorption  coefficient  spectra,  and  regressed  on  Rr5(2)  using 
the  multiple  linear  regression  described  in  equation  (2), 

ln[S275— 295]  =a  +  /?ln[jRrs(443)]  +  y  ln[Rrs(488)]  + 

(2) 

dln[Prs(53l)]  +eln[Prs(555)]  +Cln[Prs(667)] 

where  a  =  -3.4567,  £  =  0.4299,  y  =  0.0924,  3  =  -1.2649,  £  =  0.8885,  and 
£  =  —0.1025  are  the  derived  regression  coefficients.  This  data  set  covers  abroad  range 
of  water  types  found  in  river-influenced  ocean  margins,  including  high-CDOM  and 
sediment-laden  coastal  environments,  productive  waters  at  intermediate  salinities 
and  oligotrophic  waters  at  higher  salinities.  The  inclusion  of  data  from  the  Chukchi 
Sea  in  summer  (ICESCAPE)  further  demonstrates  that  the  algorithm  performs  very 
well  in  areas  where  tDOM  is  relatively  low  and  biogenic  particles  are  abundant  and 
variable.  Similar  accuracy  of  S275-295  retrievals  was  obtained  using  data  collected 
seasonally  in  the  northern  Gulf  of  Mexico  in  various  water  types  ranging  from 
riverine  to  oligotrophic,  thereby  indicating  that  the  approach  used  in  the  algorithm  is 
robust. 

MODIS  Aqua  data  were  accessed  from  the  NASA  ocean  color  website  (http:// 
oceancolor.gsfc.nasa.gov/).  Monthly  NAO  and  AO  indices  were  acquired  from  the 
NOAA  Earth  System  Research  Laboratory  website  (http://www.esrl.noaa.gov/psd/ 
data/climateindices).  The  Mackenzie  River  discharge  data  were  obtained  from  the 
Environment  Canada  Water  Office  website  (http://www.wateroffice.ec.gc.ca). 
Eurasian  river  discharge  was  obtained  from  the  ArcticRIMS  website  (http:// 
rims.unh.edu/data.shtml). 
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